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Constraining the Lorentz Factor of a Relativistic Source from its 

Thermal Emission 
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ABSTRACT 

We propose a direct and simple method to measure the Lorentz factor of a 
relativistically expanding object such as gamma-ray bursts. Only three measure¬ 
ments, i.e. the thermal component of the emission, the distance and the variable 
time scale of the light curve, are used. When the uncertainties are considered, 
we will obtain a lower limit of the Lorentz factor instead. We apply this method 
to GRB 090618 and get a lower limit of 22 for the Lorentz factor. This method 
can be used for any relativistically moving object, such as gamma-ray bursts and 
soft gamma-ray repeaters. 

Subject headings: gamma-ray: bursts, Lorentz factors 


1. Introduction 


For a relativistically moving object it is crucial to know the Lorentz factor (r). In 
the context of astrophysical objects, there are several phenomena which are claim ed t o 
be associated with relativisti c motions, such as ga mma-ray bursts (GRBs) (IPiranl 120041 1. 
active galactic nuc lei (AGN) (IGhisellini et al.l Il993t ) . and soft gamma-ray repeaters (SGRs) 
ONor ri s et ah Il99llh Especially for GRBs, the Lorentz factors of jets are believed to be about 
a few hundred. 


The identification of high T is from the compactness problem (Rees 19671 ). There are 
several methods to const rain the Lorentz factor_of the relativistic motion, like superlumi- 
nal motion for the AGN f Verm eulen & Cohen 1994), the afterglow deceleration time of the 


GRBs ( Meszaros & Reesl 1997f). the thermal emis sion from a standard fireball at the photo¬ 
sphere radius ( Nakar et al. 2005 : Pe’er et ahll2007 L the compactness problem for the observed 
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high energy phot ons (IFenimore et al.lll993l: iLithwick fc Saril 120011) . and the quiet period of 


Meszaros & Rees 

1997: 

Nakar et al. 

be obtained ( 

Fenimore et al. 

1993; 


Our method can directly derive (or constrain the lower limit of) T if the thermal emis¬ 
sion is identified, while this therm al component ha s been observed in many cases, such 
as GRB 090618 flPaee et allhoilh. G RB 090902B rtPe’er et al.l hoij l and GRB 100724B 


flGuiriec et al. 20111). 


Fan et al 


(12012 1 suggested the thermal component may even domi¬ 


nate the prompt emission in some GRBs. We mainly concentrate on the GRBs, because the 
GRB jets are mostly believed to be pointing toward the observers. Considering the beaming 
effect of the relativistic jet, spherical geometry can be safely taken only if the jet opening 


angle 9j > -L The general scenario of the GRB thermal emission is described by T oma et al. 


120 Ilf): multiple conical shells are expanding, which produce the thermal pulses, while the 
internal shocks from the collision of these shells produce the non-thermal emission at larger 
radii. The superposition of these two kinds of emissions is the observed spectrum, ffere 
we concentrate on the thermal component to derive the Lorentz factor of these shells. We 
describe the basic concept in section §2, apply the method to GRBs in section §3, and discuss 
the results in section §4. 


2. Basic concepts 

For a relativistically expanding shell, the thermal photons are coupled with the electrons 
at the beginning. With the expansion of the relativistic flow (either spherical or conical in 
shape), it becomes optically thin at the photosphere, where the optical depth drops down 
to unity. At this radius almost all the thermal photons escape, their luminosity following 
the blackbody spectrum. Beyond this radius, as no more photons are produced, emission 
drops down quickly and can be neglected, the energy stored in the plasma will be released 
by the internal or external shocks at larger radii, where nonthermal emission (synchrotron 
radiation for instance) dominates. Therefore, the observed thermal emission can be taken 
as blackbody emission at the photosphere. 0 


1 There is another possibility, namely that in the optically thick stage, photons coupled with electrons 
are no longer produce d. With the expa nsion of the relativistic shells, the shape of the spectrum becomes 
a Compton spectrum (|Vurm et al. 2013j). The number density of the photons is then smaller than from a 
blackbody emission. Consequently, the luminosity at the photosphere is lower than from a blackbody with 
the same temperature. This leads to a lower limit for the Lorentz factor. Though, in principle, this kind of 
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Because of the beaming effect, for a blackbody sphere expanding with Lorentz factor 
T (r 1), the received flux is mainly emitted from the region within 6 < yr, where 9 is 
the angle starting from the line of sight, as shown in Figure 1. The blackbody flux we can 
receive satisfies 

ff[(l + z)T oiI }H7r(^l 


™r> s < f ^ 

1 thermal — 1 bb 


4v tDI 


where -^thermal i s the observed flux of the thermal component, Fbb is the observed flux by as¬ 
suming the radiation is from a blackbody in thermal equilibrium, cr = 5.67x 10” 5 erg cm -2 s -1 1 
2 is the redshift, T 0 b s is the observed temperature, R p h is the photosphere radius of the ex¬ 
panding fireball, and D L is luminosity distance. Thermal — -^bb comes from the fact that the 
Thom pson scatte ring; photosphere can be larger than the absorption/radiation photosphere 


see 


Vurm et al. 2013, for more details). Notice that the equation is also applicable for a 


conical jet with opening angle 9j > ^ when the line of sight is inside the jet angle. 

As almost all the thermal photons in the plasma will be radiated at the photosphere 
suddenly, th e o bserved duration of this emission is generally dominated by the angular time 
scale (Piran 1999 ). 5t ang ~ where c is the speed of light. In some cases the angular 

time scale is just a lower limit, e.g., for the thick shell case before the spreading radius, the 
observed time scale is dominated by the thickness of the shell, or several thermal pulses are 
too close to be distinguished, both cases imply 5t e > 5t 


ang> 


i.e. 


5t(fs > 


(1 + z)R p h 

ZTR: ' 


( 2 ) 


Combining eqs. (H]) and (J2]), we get the lower limit for the bulk Lorentz factor: 


T > 


Fobs n2 
L thermal^L 


l(l + z)VZ2b.(2<*» 


(3) 


As all the variables on the right hand side are observables, the T is then a directly measurable 
quantity, by measuring the spectrum of the thermal component and the redshift of the object. 
The equality in the equation above applies in the ideal condition: the emitting region is 
geometrically thin, and is in thermal equilibrium. 

There is another constraint from the fact that the thermal component is mainly released 
at the radius where the shell becomes optically thin, i.e., the optical depth r = 1. This 


thermal emission can be distinguished from the blackbody emission (Plank spectrum) by the spectral shape. 
It is not easy to distinguish in real data. The combination leads to a lower limit rather than a certain value, 
which is shown in eq. ©■ Whether or not the Compton spectrum is being distinguished from the Plank 
spectrum, eq. ([!]) splits into either an equality or a lower limit. 
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determines the radius of the photosphere (IMeszaros et al.ll2002l; iDaisne &; Mochkovitchll2002f) 

(4) 


La T 

ph 8nT 3 mpC 3 


where L is the kinetic luminosity of the ejecta, which is taken as cold when the thermal 
pho tons are emitted, ctt is the Thompson cross-section and m p is the mass of the proton. 
Pe’er et al.1 (120071 ) considered this relation together with the blackbody radiation luminosity 
(i.e., eq. (JT]) ) to derive the initial Lorentz factor without considering the time scale (i.e., 
eq. ((2]))EI One should consider the consistency between the three equations: eqs. (JT]), (J2]) 
and (141) , but keeping in mind the uncertainties of eq. (JTJ) arising from the facts: 1. the 
efficiency of conversion of the kinetic energy to the observed y-ray energy is not known, it 
may include the Poynting flux making the efficiency even more uncertain; 2. the amount 
of electron-positron pairs is not known, which will affect the optical depth; 3. the matter 
should be fully ionized, otherwise another parameter characterizing the ionization should be 
introduced; 4. Eq. (|I]) is only suitable when the Thompson scattering is the dominant cause 
of opacity; 5. F t ^ s rmal = F hh is assumed. In the optimistic condition, one can combine eq. 
(Ill) and eq. (Jl]) to get the Lorentz factor by introducing an efficiency parameter Y = L/L 1 
(L 7 is the observed y-ray luminosity) (IPe’er et al.l 20071 ) 


r = 


aa T (l + zfYL.,T^ 

8n7V C 3 Dt P obs1 / 2 
OTUllpC -tVL-r thermal 


1/4 


(5) 


3. Application to GRBs 


Taking the typical values of GRBs into eq. (J3]), one gets the scaling law: 


thermal,—281 


r > 110-FL^JL, _rT)t j .28(—-—) 1 T l 2 


obs,8^©,-2- 


( 6 ) 


The conventional notation Q = Q x x 10® is used throughout this letter. The greatest 
uncertainty here is the choice of time scale ~ 0.01 s, as the variability o f the thermal 
component has not been identified in so small time bins flPage et al.l 20 lit Guiriec et al. 


2 Though they did not consider the angular time scale explicitly, when they argue that the thermal emission 

after the peak at around 8 s of GRB 970828 comes from the off-axis emission, the angular time scale shows 
the off-axis emission could not last a time of more thau_10 s. In addition, if the followed thermal component 
was coming from the high latitude as Pe’er et al. (12007 ) claimed, the temperature decay with time should 
be t~ l but not f -0 51 as observed. Therefore, we should conclude the decaying time scale of seconds is not 
from the off-axis emission. 
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20131; IBurgess et al.l 120141). Though the o bserved typical time scale for the GRB pulses is of 
the order of 0.1 s fjNakar & Piran 2002), they are believed to be from the internal shocks 
which occur farther out than the photospheric radius. Consequently, the time scale of the 
photosphere radius should be smaller. It may also be supp orted by the finding that the GRB 
pulses also contain millisecond variabilities (jWalker et al.l 120001) . 


There are several GRBs which were reported t o have the the rmal component in the 
prompt e mission: GRB 060 218 ( Gam pan a e; al. 20061). GRB 09 0618 ( Page et al. 2011), GRB 
0909 02B f Pe’er et al. 20121). GRB 100316D (Starling et alihoil ). GRB 100724B f Guiriec et al 


2 011 ). GRB 120323A (IGuiriec et al.l 2013|) and more in a recent analysis (IBurgess et al. 


2014). However, most of them have no light curve of the thermal component, which pre¬ 
vents us from estimating the time scale. Among them, GRB 120323A is a strong short burst. 
It makes the tim e reso lv ed sp ectral analysis possible, which shows a thermal time scale less 
than 0.1s fjGuiriec et al. 20131 ). Unfortunately, the redshift of this burst is not known. 


We take GRB 090618 as an example. It has been reported to have a thermal component, 


with P^rmai — 10 6 erg cm 2 s 1 , z — 0.54, D L = 9.6 x 10 27 cm, kT ohs ~ 20 keV flPage et al. 


2 01 lh . The most uncertainty is in the time scale of the thermal pulses, while T is sensitive 
to it. We choose 5t ~ 0.01s as a typical time scale for the thermal pulses. Then the Lorentz 
factor is greater than ~ 22 by taking the values into eq. (J6j) . The uncertainty of the time scale 
of the thermal component makes the constraint trivial. This constraint will be enhanced if 
the short time resolved spectral analysis can be performed, just like GRB 120323A. 

We may consider the consistency of this method with other methods in the Lorentz 
factor constraint. Taking the observed quantities into eq. (I5jl . we get the Lorentz factor to 
be 130U 1 / 4 . This is consistent with the lower limit above (keep in mind, even if eq. ([5]) is 
not consistent with other methods, that might be because eq. ([5]) is not suitable for the 
r easons lis ted abo ve). One can also estimate the Lorentz factor by the afterglow light curve 
(Meszaros & Rees 1997). From the X-ray afterglow of GRB 090618 no peak was observed, 
which indicates the peak time (deceleration time) should be smaller than 300 s, as shown in 
the light curve (IDado fe Par 2010). This upper limit of peak time (300 s) indi cates a lower 
limit for the initial Lorentz factor of roughly 150 obtained by using eq. (1) ofiLii et al. (2012) 
taking observed quantities and assuming an ambient number density 1 cm -3 . Hence the 
direct Lorentz factor is also consistent with the estimation from the aftergl ow dec eleration 
time scale. It should be consistent with the lower limit given by Lithwi ck fe Sari (12001 ). 
However, the required high energy photons from GRB 090426 were not reported. On the 
other hand, a lower limit is always consistent with another lower limit obtained by a different 
method. 


Besides the prompt emission, one can apply our method to estimate the Lorentz factor 
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from the thermal emission of the afterglow. Consequently, the Lorentz factor characterizes 
the afterglow rather than the prompt emission region. GRB 090618 has also been re porte d 


with a thermal component from the afterglo w at hundre ds o f seconds (Page et al 


2011 ). 


However, it was in the steep decline stage f lZhang et al.l 2006), which is mainly taken as 
indicating high latitude emission from the prompt radiating shell. Hence, this thermal 
component cannot be directly used to derive the Lorentz factor. 


GRBs are widely known to be highly relativistic with T ~ 100 (lPiranl 2004J). For a given 
T ~ 100, we can see from eq. (EJ), that the temperature T 0bs must be relatively low (say 
roughly < lOkeV), otherwise, the thermal radiation will be extremely high and cannot be 
missed by the spacecraft. The relatively low temperature might be the reason why so few 
GRBs have been firmly determined to have the thermal component. 

Notice that very few GRBs prompt emission are identified with a thermal component, 
and the thermal components are all identified with observed temperatures less than 10 keV. 
The main reason is that the current on-boading satellite missions mainly cover relatively 
high energy bands (more than tens of keV). It may b e resolved by the lower band soft j - 


ray teles copes scheduled to be launched, like SVOM (Gotz fe SVOM Collaboration! 120121 1 . 


JANUS (Falcone et al. 2009) and LOFT (Amati et al. 


2015). After identifying the thermal 


component, one needs to identify the time scales of the thermal pulses, which means the 
time bin should be small enough to perform a spectral analysis. This is still challenging but 
might be available for strong soft bursts in the near future. Otherwise, only lower limits can 
be achieved. 


Discussion 


To clarify, we emphasize the scenario for the thermal and non-thermal components in 
the following. The relativistic warm shells produced from the central engine are ejected 
stochastically. With the shells expanding spherically the optical depth for the trapped ther¬ 
mal photons decreases. At the transparent radius, the shell emits almost all the thermal 
photons, which are observed as the thermal component. Each shell will produce one thermal 
pulse with a short time scale. The combination of any two shells produces one internal shock, 
which generally occurs at a larger radius, and produces the non-thermal component. This 
non-thermal component has a longer time scale. The superposition of these two components 
is what we observe. For the purpose of obtaining the Lorentz factor, we need to extract the 
thermal component. For GRBs, the thermal component should never be dominant. Oth¬ 


erwise, there will be more energy le: 
problem and the efficiency problem l 


tAn the plasma. This is related to the energy budget 


Piran 2004j). This makes the extraction of the thermal 
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light curve more difficult. 


In general, this method can be applied to any object radiating thermal emission, and 
with a relativistic velocity. There are several other objects to which this method may be 
potentially applicable, though shortcoming s are obvious at the present stage. SGRs are also 
thought to be relativistic events flCheng fc Wang| |2003l. for example). However, the bl ackbody 


emission in this case is mainly taken to be generated from the progenitor flHurlev et al. 


2005 ). This method cannot be directly applied to AGN, as AGN jets are mainly pointing 
to ot her directi ons, they are all continuous jets, and there is no confirmed thermal radiation 
(IFossati e_t al. 119981) . X-ray binaries are believed to have relativistic jets too. However, the 
number of on-axis case is very limited, and the thermal component is not easily distinguished 
as arising from the jet or from the disk. Tidal disruption events also have mildly relativistic 
ejecta. But no associated thermal component has been identified. The low Lorentz factor 
may be beyond the restriction 1/T < 9j, and is not easy to obtain by our method. This 
method may also be applied to the GRB afterglow stage if there is a thermal component. 


We should also notice that it is applicable only if the following conditions are satisfied: 

1. The jet opening angle should be greater than 1/T, otherwise the shortest time scale is not 
R p h/2T 2 c. Fortunately, for the GRBs the jet opening angle is generally greater than 1/T. 

2. The line of sight should be inside the jet opening angle, and this is always the case for 
GRBs. 3. Though in the GRB fireball scenario, it is unlikely the ejecta is still accelerating 
at the photosphere radius, it is possible in general. Then the derived Lorentz factor is a 
lower limit again. 


One can also constrain the Lorentz factor by other methods, especially the upper limits. 
By combining several methods, once the Lorentz factor is constrained in a very narrow range 
without inconsistency, that means the equality of eq. (EJ) can be applied. This indicates that 
for such a burst the thermal emission comes from a thin shell and no heating of the electrons 
occurred. 
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Black-body radiating shell 


Fig. 1.— Sketch for the conical radiating shell at the photosphere, which is moving with 
Lorentz factor T. Because of the beaming effect, roughly only the region inside 1/T can be 
seen by the observer. This determines the time scale of the thermal pulse. 



